Physical properties of FeSeo.sTeo.s single crystals grown under different conditions 
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We report on structural, magnetic, conductivity, and thermodynamic studies of FeSeo.sTeo.s single 
crystals grown by self-flux and Bridgman methods. The lowest values of the susceptibility in the 
normal state, the highest transition temperature T c of 14.5 K, and the largest heat-capacity anomaly 
at T c were obtained for pure (oxygen-free) samples. The critical current density j c of 8.6 x 10 4 A/cm 2 
(at 2 K) achieved in pure samples is attributed to intrinsic inhomogeneity due to disorder at the 
anion sites. The samples containing spurious phase of FesCU show increased j c up to 2.3 x 10 5 A/cm 2 
due to additional pinning centers. The upper critical field H C 2 of ~ 500 kOe is estimated from the 
resistivity study in magnetic fields parallel to the c-axis. The anisotropy of the upper critical field 
7if C 2 ~ H ra ^/H c 2 reaches a value ~ 6 at T — > T c . Extremely low values of the residual Sommerfeld 
coefficient for pure samples indicate a high volume fraction of the superconducting phase (up to 97 
%). The electronic contribution to the specific heat in the superconducting state is well described 
within a single-band BCS model with a temperature dependent gap Ao = 27(1) K. A broad cusp- 
like anomaly in the electronic specific heat of samples with suppressed bulk superconductivity is 
ascribed to a splitting of the ground state of the Fe 2+ ions at 2c sites. This contribution is fully 
suppressed in the ordered state in samples with bulk superconductivity. 



INTRODUCTION 



The recent finding of high-temperature superconduc- 
tivity in iron pnictides [H-Q and iron chalcogenides Q 
inspired an immense research activity in complex mate- 
rials similar to the discovery of high-T c cuprate super- 
conductors two decades ago. Among the new iron-based 
superconductors, the iron chalcogenides are of particular 
interest due to a simple crystal structure (Fig. 1). It con- 
sists of Fe ions tetrahedrally coordinated by Se and Te 
arranged in layers stacked along the c-axis in the tetrag- 
onal lattice without any other interlayer cations, as for 
example in pnictides. Therefore, the iron chalcogenides 
were believed to be most suitable for investigation of the 
interplay of competing orders and pairing mechanisms. 
However, recently it was recognized in Fe(Se,Te) that Fe 
ions can occupy also the 2c positions in the anion plane 

Nearly stoichiometric FeSe becomes superconducting 
below 8 K at ambient pressure [3] but the transition tem- 
perature can be enhanced up to 37 K by application of 
external pressure [3, l8| . The superconducting properties 
of FeSe were shown to be extremely sensitive to devia- 
tions from the stoichiometry 0, [l(| ■ The substitution of 
Se by Te in FeSe increases the temperature of the super- 
conducting transition to a maximum of about 14 K for 
50 % of replacement [H, 12 1. The intrinsic disorder due 
to random distribution of the Se and Te ions among the 
anion sites in the crystal lattice is expected to reduce the 
sensitivity of FeSe to non-stoichiometry. Subsequently, 
several authors recently reported on the successful prepa- 
ration of high quality single crystals of FeSei _^Te^ using 
Bridgman \1ML^ and flux methods [IE Ed]. Tlle exist " 
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FIG. 1. (color online) Crystal structure of FeSeo.sTeo.s 
demonstrating the tetrahedral coordination of Fel ions in 2a 
positions (blue spheres) by anions of Se and Te (green spheres) 
and the interstitial Fe2 ions in 2c position (red spheres) in the 
anion plane. 



ing data show large variations in magnetic and conduct- 
ing behavior of Fe-Se(Te) obviously related to differences 
in the growth conditions and purity of the starting ma- 
terials. Hence, only little information about the intrinsic 
properties of FeSeo.sTeo.s is available so far. 

Here we report on the magnetic susceptibility, resis- 
tivity and heat-capacity measurements of FeSei-^Tea; 
with nominal concentration x = 0.5. The samples were 
prepared in various ways and exhibit significantly differ- 
ent behavior which was investigated in the temperature 
range 2 - 400 K and in external magnetic fields up to 
140 kOe. 
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EXPERIMENTAL 



Single crystals of FeSeo.sTeo.s were grown by self-flux 
and Bridgman methods. Flakes and shots of high-purity 
elements, 99.99 % Fe, 99.999 % Se, and 99. 999 % Te, 
were used in the growth experiments. To get a mini- 
mal amount of oxide impurity we additionally purified 
Se (for growth runs F216) and Te (for growth runs F213 
and F216) by zone melting. Handling of these samples 
was done in an argon box with residual oxygen and water 
content less than 1 ppm. The growth of single crystals 
was performed in evacuated double quartz ampoules. For 
different runs in the self-flux method the soaking temper- 
ature varied between 900 and 1100° C. The cooling rate 
varied between 1 and 60 ° C /h. In the Bridgman method 
we used pulling rates between 0.5 and 2 mm/h and rota- 
tion speed of 2-5 turns/min. Final treatment of the sam- 
ples in both preparation methods was done at 410° C 
for 70 - 100 hours followed by quenching in ice water. 
The plate-like samples with shiny faces with dimensions 
~3x5 mm 2 in the ab plane and thickness up to 0.5 mm 
along the c axis were obtained on the top of the solidified 
boule. The single-crystallinity of the grown samples was 
checked by the single crystal X-ray diffraction. The sam- 
ple composition was investigated by Energy Dispersive 
x-ray (EDX) analysis. The phase content of the samples 
was analyzed by x-ray powder diffraction (Cu K a radia- 
tion, A = 1.540560 A) on crushed single crystals using a 
STADI-P powder diffractometer (STOE & CIE) with a 
position sensitive detector. 



The magnetic measurements were performed in a tem- 
perature range 2 - 400 K and in magnetic fields up to 
50 kOe using a SQUID magnetometer MPMS 5 (Quan- 
tum Design). Electron-spin resonance (ESR) measure- 
ments were performed in a Bruker ELEXSYS E500 CW 
spectrometer at X-band frequencies (9.36 GHz) equipped 
with a continuous He gas-flow cryostat in the tempera- 
ture region 4.2 < T < 300 K. The heat capacity was 
measured by a relaxation method using a Quantum De- 
sign physical properties measurement system (PPMS) in 
a temperature range 1.8 - 300 K and magnetic fields up 
to 90 kOe. The magnetic field was applied parallel to 
the c-axis of the samples. The resistivity studies were 
done on rectangular samples by a four-point method in 
the temperature range 2 - 300 K and in magnetic fields 
up to 140 kOe using a He-flow cryostat (Oxford Instru- 
ments) and using the resistivity measurement option of 
the PPMS in magnetic fields up to 90 kOe. The contacts 
were made by conductive silver paint. 
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FIG. 2. (color online) a) Room temperature x-ray diffrac- 
tion patterns for FeSeo.sTeo.s samples prepared under differ- 
ent conditions. Stars indicate the impurity phase; b) Exper- 
imental (open circles) and refined (black line) x-ray diffrac- 
tion patterns for sample F213. Blue line- difference between 
experimental and calculated intensities. Vertical lines mark 
Bragg positions: top row- main tetragonal phase; bottom- 
impurity hexagonal phase. Growth conditions: Sample BrN5 
(Bridgman): soaking temperature T s = 1084° C; soaking time 
t s =24 h; puling rate 0.5 mm/h. Sample F213 (self-flux): 
T S =1100°C; t s =72 h; cooling rate 60°C/h. Sample F216 
st. 1 (self-flux): T S =1100°C; t s =72 h; cooling rate 60°C/h. 
Sample F216 st. 2 (self-flux): T s =1100 ° C; t s =72 h; cooling 
rate l°C/h. 



EXPERIMENTAL RESULTS AND DISCUSSION 
X-ray analysis 

The x-ray diffraction patterns for the samples under in- 
vestigation (Fig. 2a) are consistent with tetragonal sym- 
metry P4/nmm for the main FeSeo.sTeo.s phase and with 
hexagonal symmetry P63/mmc for the impurity phase 
Fe7Se§. The refinement was performed using the FULL- 
PROF SUITE [3. The results of the refinement for one 
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of the samples are shown in Fig. 2b. The occupation of 
Te and Se at the 2c sites was refined constraining their 
sum to unity in agreement with the EDX data. A similar 
constraint was used for the occupation of Fe ions for the 
main phase allowing for two different sites (2a and 2c). 
The EDX measurements on samples from various growth 
experiments find no indication that the different methods 
used for the crystal growth yield different compositions. 
No significant gradient in the Te/Se ratio along the sam- 
ple was observed. The Se content for samples from differ- 
ent growth runs varied from 0.455(6) to 0.498(10), and 
the respective variations of the Te content were within 
the range 0.555(4) to 0.502(10). The observed maximal 
deviations of Se/Te ratio for some samples of about 10 % 
from the nominal 50/50 % composition, correspond to 
FeSeo.45Teo.55. The ratio Fe/(Se+Te) in the samples was 
deviating less than 2 % from the stoichiometric value 
indicating a low content of excess iron. The results of 
the refinement for different samples are given in Table 
1. The refined occupation factors for Te and Se ions for 
the main tetragonal phase are in general agreement with 
the EDX data. The occupation of iron on the 2c sites 
was ~ 6 — 9 %. The values of the lattice parameters for 
different samples varied within the range 3.800- 3.803 A 
for the a(b) lattice constant, and 6.026 - 6.047 A for the 
c parameter and are close to those reported in Ref. 13 
for single crystals of similar composition. 



Susceptibility and magnetic hysteresis 

Fig. 3a shows the temperature dependence of the mag- 
netic susceptibility \ for samples from different growth 
runs measured on cooling in a field of 10 kOe applied 
parallel to the c-axis. The samples prepared from the 
non-purified Se and Te and handled in air have a rel- 
atively high susceptibility, ~ (1 — 2) x 10~ 2 emu/mol. 
Below 300 K the susceptibility continuously increases on 
decreasing temperature, then shows a pronounced down- 
turn at 125 K and a clear anomaly at the superconduct- 
ing transition at T c ~ 14 K. The downturn at 125 K 
is related to the Verwey transition of Fe3C>4 I19|. w hich 
was also evidenced by electron-spin resonance [20( and 
specific-heat measurements (see below) . The ESR exper- 
iment did not reveal any intrinsic absorption of localized 
moments. 

Samples prepared from the purified elements and han- 
dled in the argon box exhibit susceptibility by one or- 
der of magnitude lower and do not show any anomaly at 
125 K. The susceptibility of the purest oxygen-free sam- 
ple, F216 step 1, prepared by fast cooling (l°C/min) 
manifests a non-monotonous temperature behavior with 
a broad maximum at around 180 K. It also shows the 
sharpest anomaly at around 14 K at the superconducting 
transition. Surprisingly, repeated melting of this sample 
followed by slow cooling with the rate of l°C/h, dras- 



tically changes the magnetic behavior. The susceptibil- 
ity of this sample is temperature independent at high 
temperatures and shows a paramagnetic tail at low tem- 
peratures which might be associated with localized iron 
ions. Recent investigations of FeSe(Te) samples with iron 
excess [1, [13, 21 [ established a localized magnetic iron 
moment at the 2c sites which can suppress superconduc- 
tivity. However, both well and poorly superconducting 
samples contain comparable amount of the Fe ions at the 
2c sites (Table 1). An explanation of the low-temperature 
susceptibility tail in the poorly superconducting samples 
due to the hexagonal-phase can be excluded because the 
amount of the hexagonal phase can be strongly reduced 
in the oxygen-free samples. The apparent effect of the 
hexagonal impurity on the susceptibility can be revealed 
by comparing the oxygen-free samples F213 and F216 
step 1, containing respectively 5.5 and 1.4 % of the hexag- 
onal phase. 

Fig. 3b shows the temperature dependence of the zero- 
field cooled (ZFC) and field-cooled (FC) susceptibility 
measured in a field of 10 Oe applied parallel to the c- 
axis. The measurements along the a&-plane showed sim- 
ilar results except that the diamagnetic response was sub- 
stantially higher for the field parallel to the c-axis. The 
samples containing both, oxide and hexagonal-phase im- 
purities manifest rather broad transitions into the super- 
conducting state. At the same time, the temperature 
of the onset of the superconductivity, T° n , shows only 
a small change at around 13.5 K for these samples. We 
found the sharpest transition and, respectively, the high- 
est T c of 13.8 K (with T° n of 14.5 K) for samples grown 
from the purified elements under fast cooling conditions 
only. The values of 4irx at 2 K calculated from the data 
shown in Fig. 3b are far above unity indicating that the 
low-field susceptibility of samples with high diamagnetic 
response is dominated by demagnetizing effects. The vol- 
ume fraction of the superconducting phase was therefore 
estimated from \z fc 011 needle-like crystals for the field 
applied along the a&-plane with the smallest demagne- 
tizing factors and reached 98 % for the purest samples. 
At the same time, in the poorly superconducting sample 
the onset of the transition remains at the same tempera- 
ture as in well superconducting samples, but it is masked 
by a temperature independent paramagnetic background 
and therefore can be evidenced only on high magnifica- 
tion of the data as shown in the inset in Fig. 3b. The 
volume fraction of the superconducting phase with T° n 
of 14.5 K in this sample is by more then two orders lower 
than in the well superconducting sample suggesting the 
filamentary superconductivity. 

Fig. 4 presents the magnetization hysteresis loops for 
different samples measured at 2 K with the field ap- 
plied along the c-axis. The hysteresis loops for samples 
with a high volume of the superconducting phase have a 
symmetric character which indicates dominant bulk pin- 
ning and small contribution of surface pinning [22j. Be- 
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TABLE I. Structural data obtained from the Rietveld refinement for FeSeo.sTeo.s samples 



Sample 


Occupation 


Lattice constant 


Tetragonal 
phase 
[%] 


Hexagonal 
phase 
[%] 






Chi 2 


Fel 
/2a/ 


Fe2 

/2c/ 


Se 
/2c/ 


Te 
/2c/ 


a, b 

[A] 


c 

[A] 


Br N5 


0.907(3) 


0.093(3) 


0.48(1) 


0.52(1) 


3.8020(4) 


6.0489(9) 


94.5 


5.5 


3.53 


2.48 


2.02 


F213 


0.930(5) 


0.070(5) 


0.49(1) 


0.51(1) 


3.8011(2) 


6.0409(7) 


94.5 


5.5 


3.68 


2.54 


2.10 


F216 step 1 


0.929(3) 


0.071(3) 


0.49(1) 


0.51(1) 


3.8025(3) 


6.0300(9) 


98.6 


1.4 


3.0 


2.36 


1.62 


F216 step 2 


0.937(3) 


0.063(3) 


0.52(1) 


0.48(1) 


3.8000(3) 


6.0257(8) 


97.9 


2.1 


3.06 


2.49 


1.51 



low the transition temperature the diamagnetic response 
in these samples dominates over the full measured field 
range. The samples containing oxide impurities show 
an enhanced (by ~ 20 %) magnetization compared to 
the oxygen-free samples due to additional pinning cen- 
ters of Fe3C>4. A strongly contrasting hysteretic behavior 
was observed for the oxygen-free samples prepared under 
slow-cooling conditions. These samples exhibit a con- 
siderably reduced width of the hysteresis loops, by two 
orders of magnitude, suggesting full suppression of the 
bulk superconductivity. 

Surprisingly, such a drastic change of the magnetic 
properties of the sample with suppressed bulk supercon- 
ductivity is not accompanied by notable changes in the 
structural data. The x-ray diffraction patterns of these 
samples show a comparable amount of the tetragonal and 
impurity phases (Fig. 2a). Their refinement data (Table 
1) do not reveal any essential variations of the site oc- 
cupancy for the iron ions and for the anions as well as 
of the lattice parameters of the main tetragonal phase. 
Although any comparison of the low-temperature super- 
conducting data and the x-ray diffraction data taken at 
room temperature is not plausible, our results indicate a 
rather subtle role of the structure which probably cannot 
be resolved using conventional x-ray facilities. 

Figures 5a and b show, respectively, the magnetic hys- 
teresis and critical current at different temperatures for 
the purest samples with bulk superconductivity. The 
critical current density was estimated from the width of 
the hysteresis loops using the Bean model for hard su- 
perconductors [U, |3] . At 2 K the critical current den- 
sity j c at zero field reaches a value of 8.6 x 10 4 A/cm 2 . 
For the sample grown by Bridgman method we obtained 
j c = 9.4x 10 4 A/cm 2 . For the sample grown from self-flux 
and intentionally oxidized after the first step we found 
nearly 3 times higher values of j c = 2.3 x 10 5 A/cm 2 . 
Thus, the oxidation allows to increase the critical cur- 
rent. The critical current for the purest sample with high 
volume fraction of the superconducting phase decreases 
approximately three times in fields up to 20 kOe but then 
flattens at this level up to the largest measured fields sug- 
gesting a high current-carrying ability of the material. 
An estimation of critical currents for T = from fits to 
the experimental data was performed using a generalized 
power-law dependence j(T) = j(0)[l - (T/T c ) p ] n , with 



p = 0.5, n = 1.5 and T c — 13.8 K. The respective fit is 
shown by the dotted line in the inset of Fig. 5b, yielding 
a value j(0) = 1.7 x 10 5 A/cm 2 . The calculated value 
of j(0) is close to those determined for high-quality su- 
perconducting single crystals of Ba(Fei_ a; Co a ;)2As2 [25j ]. 
Table 2 summarizes the critical current densities calcu- 
lated from the hysteresis loops at 2 K together with the 
critical temperature T c and lower critical field H c \ deter- 
mined from the magnetic data. 



Electrical resistivity 

Fig. 6a presents the temperature dependences of the 
resistivity for the studied samples. The oxygen-free sam- 
ples have a lower value of the resistivity in the normal 
state compared to the samples with oxide impurity which 
can be attributed to increased scattering of charge carri- 
ers on additional impurity centers related to Fe304. The 
resistivity of samples with high volume fraction of the 
superconducting phase exhibits a metal-like temperature 
dependence below 200 K down to T c as shown in the in- 
set of this figure. Such a behavior was established earlier 
only for samples with a low amount of excess iron 2l| . In 
contrast, the resistivity of samples with suppressed bulk 
superconductivity shows a temperature independent be- 
havior at high temperatures down to T c . The resistivity 
of these samples drops at approximately the same tem- 
perature as in the samples with high superconducting 
parameters, but the resistive transition is strongly broad- 
ened and shifted to lower temperatures in agreement with 
the susceptibility data. 

Fig. 6b illustrates the effect of magnetic field on the 
in-plane resistivity in the transition region for one of the 
oxygen- free samples (F213). The field was applied par- 
allel and perpendicular to the c-axis, and the measure- 
ments were done on warming after cooling in zero field. 
The resistivity curves are displaced to lower tempera- 
tures with increasing magnetic field with a stronger shift 
of the transition for the field parallel to the c-axis than 
for the perpendicular configuration. The obtained data 
are in general agreement with the respective results on 
the superconducting FeSei-^Tea, samples with different 
x reported earlier 17, 26| . 

In Fig. 7 the temperature dependences of the upper 
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TABLE II. Superconducting parameters: transition temperature, critical current, upper and lower critical fields for FeSeo.sTeo.s 
samples 



Sample 


rpon 

■*• c 




if 










(K) 


(kA/cm 2 ) 


(kA/cm 2 ) 


(kOe) 


(kOe) 


(kOe) 






(2 K) 


(0 K) 


(2 K) 


[H || c] 


[H || ab] 


Br N5 


13.5 


94 


250 


3 


510/1300* 


850 


F213 


14.5 


77 




1.4 


570/1300* 




F216 step 1 


14.5 


86 


170 


1.2 


490/1500* 




F216 step 2 


7.6 


0.8 




0.1 







* Estimated from the heat capacity measurements. 



critical field H C 2 (T) estimated using the criterion of 50 % 
drop of the normal state resistivity p n are shown. The 
calculated data show a similar behavior for different sam- 
ples with the shift on the temperature scale correspond- 
ing to a difference in their transition temperatures. On 
approaching T c the slope of the H C 2(T) curve for the 
configuration H \\ c becomes smaller compared to that 
of the lower temperatures. Contrary, the slope of H c2 (T) 
for the configuration H _L c shows a slight increase on 
approaching T c . The anisotropy of the upper critical 
field defined as jh c2 = H&fHdZ shows a notable increase 
on approaching T c from jh c2 = 2-15 at T/T c — 0.91 to 
jh c2 = 3-6 at T/T c = 0.968 and finally reaches jh c2 = 6 
at f/T c = 0.996 (see inset of Fig. 7). This observation is 
in disagreement with the nearly isotropic behavior of the 
upper critical field reported by other studies on supercon- 
ducting FeSe(S)i_ :E Te :r with different levels of substitu- 
tion (l~71 . [26| . The higher anisotropy value of 7# c2 may re- 
flect the better quality of our samples compared to earlier 
reports. An estimation of the upper critical field H C 2(0) 
for T = K was achieved within the Werthamer-Helfand- 
Hohenberg (WHH) model for conventional superconduc- 
tors, defined by H c2 (0) = -0.69T c (dH c2 (T)/dT) \ Te for 
a weak-coupling regime [27j. The calculated data are 
presented in Table 2. The values of noH^i)) are rather 
similar for all samples, are independent of the impurity 
content, and vary in the range 490 - 570 kOe for config- 
uration of the magnetic field parallel to the c-axis. The 
estimated value of poH C 2{Q) for one of the measured sam- 
ples along the ab plane equals ~ 850 kOe. Such high 
values of H C 2(0), which are far above the Pauli paramag- 
netic limit hqH p = 1.84T C ~ 260 kOe can be attributed 
to an enhanced impurity scattering from the Fe at 2c 
sites as proposed in Rcf. 14. 



Specific heat 

Fig. 8 shows the temperature dependences of the 
specific heat for samples with different superconducting 
properties. They exhibit a similar behavior and similar 
values of the specific heat at temperatures above 15 K. 
A small anomaly at 125 K is observed for sample Br N5 



which contains Fe304. By scaling the entropy involved 
in this anomaly with that of the specific-heat anomaly at 
the Verwey transition measured in Fe3C>4 we estimated 
the amount of the oxide impurity in this sample to be 
at a level of ~ 1 mol%. The inset in Fig. 8 illustrates 
the specific heat in the representation C/T vs. T in the 
low temperature range. The samples with bulk supercon- 
ductivity exhibit a pronounced anomaly at around 14 K. 
As the transition in the poorly superconducting sample 
is smeared out, no discernible anomaly in the specific 
heat could be detected. The data in Fig. 8 indicate that 
the superconducting contribution to the specific heat is 
small compared to that of the lattice contribution which 
dominates the total specific heat. 

Fig. 9 shows the temperature dependences of the spe- 
cific heat plotted as C/T vs. T 2 for different samples are 
shown at temperatures below 5 K. A fit to the experimen- 
tal data in the range below 4.5 K by using C/T = j+f3T 2 
allows to estimate the values of the Sommerfeld coeffi- 
cient 7 related to the electronic contribution, and the 
prefactor /3 which characterizes the lattice contribution 
to the specific heat in a simple Debye approximation. For 
samples with non-superconducting behavior this proce- 
dure gives an estimate of 7 in the normal state, 7„, while 
for superconducting samples it yields the residual j r . The 
calculated values of these parameters are given in Ta- 
ble For the oxygen- free bulky superconducting sam- 
ples we obtained 7 r = 0.82 - 0.96 mJ/mol K 2 . The ob- 
tained values of the residual j r for the oxygen-free super- 
conducting samples arc mud) lower than report<'d_thus 



far by other authors for similar compositions [13l 121 



These extremely low values of "f r confirm the high pu- 
rity of our oxygen-free samples. The obtained values of 
the residual "f r for pure samples indicate that the vol- 
ume fraction of the superconducting phase reaches ~ 95 
- 96 % which agrees with the estimate obtained from the 
susceptibility data. A larger value of j r — 5.2 mJ/mol K 2 
for sample Br N5 can be probably attributed to magnetic 
contribution from the Fe304 impurity although due to in- 
sulating nature of Fe304 at low temperatures one should 
not expect any contribution to j r . However, one cannot 
exclude, for example, the presence of glassy-like disorder 
in Fe3C"4 with low energy excitations linear in temper- 
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T(K) 

FIG. 3. (color online) a) Temperature dependence of the field- 
cooled susceptibility of FeSeo.sTeo.s measured in a field of 
fO kOe applied along the c-axis. The arrow at 125 K marks 
the anomaly related to Fe3C>4 ; b) Temperature dependences 
of ZFC and FC susceptibilities of samples with bulk super- 
conductivity measured in a field of 10 Oe applied along the 
c-axis. Arrow indicates the temperature of the onset of the 
superconducting transition T° n — 14.5 K for the purest sam- 
ples. Inset: ZFC and FC susceptibilities of the sample with 
filamentary superconductivity shown on enlarged scale. 

ature or additional magnetic excitations yielding a T 3 / 2 
dependence of the specific heat. This problem needs ad- 
ditional study. For the samples with suppressed bulk su- 
perconductivity we obtained j r = 19.3 mJ/mol K 2 . We 
must also note that the calculated values of the prefactor 
j3 for samples with high superconducting parameters and 
for those with reduced superconductivity are very close to 
each other. This indicates that the electronic supercon- 
ducting contribution to the specific heat in the range of 
temperatures used for fitting has only a minor influence 
and does not much affect the accuracy of calculations. 
This justifies the above-mentioned fitting procedure and 
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FIG. 4. (color online) Hysteresis loops measured at 2 K with 
the magnetic field applied along the c-axis for different sin- 
gle crystalline FeSeo.sTeo.s samples. The data for the slowly 
cooled sample are magnified by a factor of 20 for clarity. 

the correctness of the obtained parameters. 

To get an additional independent estimate for the elec- 
tronic and lattice contributions to the specific heat we 
used the following approach. For a description of the 
phonon spectrum a combined Einstein-Debye model was 
employed. The tetragonal unit cell of FeSe(Te) with 
space group P4/nmm contains two formula units giving 
rise to a total of 12 normal modes of vibrations. Their 
contribution was simulated by two Debye terms Co and 
one Einstein term Ce with equal distribution of the spec- 
tral weight between the Debye and Einstein terms. These 
assumptions are in rough agreement with the results of 
the experimental study of the phonon density of states 
by nuclear inelastic scattering [29j and neutron scattering 
[301 ] on related Fei +a: Se superconductors. 

The characteristic Debye and Einstein temperatures, 
Go and Oe, were the input parameters for a fit to the 
experimental temperature dependence of the total spe- 
cific heat above T c by the expression 

c = c m (e m ) + c D2 {e D2 ) + c E {Q E ) + it 

The fitting parameters were varied till the minimal de- 
viations from a constant j n value in a maximal temper- 
ature range (up to 200 K) for the superconducting sam- 
ples were achieved. The temperature dependence of the 
simulated lattice specific heat with these optimized val- 
ues of 9 m = 127 K, Q D2 = 235 K and e E = 315 K 
is shown by the dotted line in Fig. 8. For the sample 
with suppressed bulk superconductivity by the dashed 
line Fig. 8 also presents a curve calculated by subtract- 
ing the normal state electronic contribution 7„T (with 
7„ = 23 mJ/mol K 2 ) from the measured total specific 
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FIG. 5. (color online) a) Hysteresis loops at different tem- 
peratures measured for the field applied along the c-axis for 
FeSeo.sTeo.s sample F216 step 1 grown by self-flux method; 
b) Critical current density j c vs. magnetic field at different 
temperatures for the same sample. The inset shows the tem- 
perature dependence of the critical current at zero field with 
the fit by the power law (dashed line). 
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FIG. 6. (color online) a) Temperature dependences of the re- 
sistivity for different FeSeo.sTeo.s samples measured on cool- 
ing in zero external magnetic field. Inset: Temperature de- 
pendences of the resistivity for samples BrN5 and F2f6 step 
I measured up to 200 K; b) Temperature dependences of the 
in-plane resistivity for oxygen-free sample F2I3 measured at 
various magnetic fields applied parallel to the c-axis. The 
arrow shows the direction of increasing field. 



heat. At temperatures above 30 K both curves nicely 
coincide with deviations less than 2 % in the complete 
temperature range up to 300 K. This consistency of the 
data again justifies the model used for the simulation of 
the phononic contribution. Note that no scaling of the 
phononic contribution for samples with high and low su- 
perconducting parameters was necessary. We found that 
the estimated values of the Sommerfeld coefficient in the 
normal state j n vary in the range 23 - 26 mJ/mol K 2 for 
different samples and are much lower than those reported 
previously for FeSei-^Tea; by other authors [l3l [ill. |28||. 
To our opinion, the reason of this discrepancy is related 
to different estimates of the lattice contribution (3ll |. 
The electronic specific heat C e for all samples was de- 



termined by subtracting the calculated lattice contribu- 
tion from the total measured specific heat. The depen- 
dences of the electronic specific heat in the representation 
C e /T vs. T are shown in Fig. 10 for different samples in 
a temperature range around the transition temperature. 
All samples with bulk superconductivity exhibit a sharp 
anomaly in C e at T c . The magnitude of the anomaly 
at T c correlates with the amount of the superconducting 
phase in these samples. In the sample with suppressed 
superconducting properties no sharp anomaly at T c is 
evidenced, but instead a broad cusp in C e centered at 
around 10 K develops. A magnetic field of 90 kOe fully 
suppresses the residual filamentary superconductivity in 
this sample, as show the susceptibility measurements, but 
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FIG. 7. (color online) Temperature dependences of the up- 
per critical field H C 2 determined using the criterion of 0.5R n 
for samples grown by different methods with magnetic field 
H parallel to the c-axis (closed symbols) and parallel to the 
ab plane (open triangles). The open squares represent the 
data for sample F213 calculated from the specific heat (see 
text). The inset shows the temperature dependence of the 
anisotropy of the upper critical field jb c2 = /H^2 i n the 
vicinity of T c for sample Br N5. 



has a negligible effect on the specific heat. This suggests 
that the broad cusp in the electronic specific heat is not 
related to the superconducting behavior. Importantly, 
a part of this cusp-like anomaly is also evidenced in C e 
for samples with bulk superconductivity as a broad right 
wing at temperatures above T c . To check whether the 
broad cusp anomaly can be an artifact due to a possible 
incorrect description of the lattice specific heat at low 
temperatures, the contribution of the lowest phonon was 
varied up to the limiting value when it coincides with 
the total specific heat. However, even with this overes- 
timated lattice contribution, the broad cusp in the elec- 
tronic specific heat still remains present for all samples, 
although with somewhat reduced amplitude (by ~ 20 %) 
and with slight shift (by ~ 1 K) to lower temperatures. 
These results indicate that the broad anomaly in the elec- 
tronic specific heat is not an artifact due to an inade- 
quate modeling of the phonon density of states, but is a 
feature reflecting the intrinsic properties of the studied 
samples. The typical Schottky-like appearance of this 
anomaly suggests an electronic origin, and the indepen- 
dence of magnetic field indicates its relation to the orbital 
degree of freedom. We speculate that it could originate 
from a splitting of the ground state of Fe 2+ ions either 
by a crystal field or due to spin-orbital coupling. There- 
fore this anomaly was simulated within a simple model 
of a two-level system. The results of the calculation are 
shown by a solid line for the sample with suppressed su- 
perconductivity using the data for a field of 90 kOe. We 



FIG. 8. (color online) Temperature dependences of the 
specific heat for different samples. The star indicates the 
anomaly at the Verwey transition due to presence of Fe304 
impurity in sample Br N5. The dashed and dotted lines 
present the lattice specific heat as described in the text. The 
inset shows the specific heat in the representation C/T vs. T 
in the transition region. 
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FIG. 9. (color online) Temperature dependences of the spe- 
cific heat in the representation C/T vs. T 2 for different sam- 
ples. 



arrived at a reasonable description of the cusp at tem- 
peratures above 7 K with the value of the ground state 
splitting Acf — 24 cm -1 and the amount of magnetic 
species corresponding to ~ 7 % per mole. The only value 
that correlates with this quantity in our samples is the 
concentration of Fe ions on the 2c sites as resulting from 
the refinement of the x-ray data. As already mentioned 
previously, these 2c site Fe ions can have a local magnetic 
moment. Therefore it is natural to assume that they can 
also be responsible for the broad cusp-like anomaly in 
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the electronic specific heat. However, the presence of this 
feature in the electronic specific heat in the samples with 
bulk superconductivity just above T c indicates that the 
magnetic moment of these Fe ions is not the main factor 
that destroys superconductivity. A remarkably sharp be- 
havior of the specific heat just below T c in samples with 
bulk superconductivity along with extremely low values 
of their residual 7 r indicates that this contribution is fully 
suppressed in the superconducting state. 




F213 

F2 16 step 1 
F216step2 
Br N5 
Schottky 
BCS F213 
BCS F216 
BCS BrN5 



7(K) 

FIG. 10. (color online) Temperature dependences of the elec- 
tronic specific heat in representation C e /T for different sam- 
ples. The solid lines represent the fits describing, respectively, 
the superconducting specific heat for the bulk superconduct- 
ing samples within the single-gap BCS model as given in the 
text, and the Schottky anomaly in the non-superconducting 
sample (F216 step 2). 

The temperature dependence of the electronic super- 
conducting specific heat was analyzed within a BCS de- 
rived a-model [3J, |35[ with a temperature dependent su- 
perconducting gap A, using a similar approach as de- 
scribed in Ref. 36 used for analysis of the specific heat in 
related Ba(K)Fe2As2 pnictides. In Fig. 10 fitting curves 
are shown by solid lines for three samples with high su- 
perconducting parameters. The fitting curves reasonably 
describe the superconducting specific heat, except in the 
range below 5 K which can be related to effects of resid- 
ual impurities. The values of the superconducting gap 
at K, Ao, derived from the analysis, vary within the 
range 26 - 28 K for different samples. The value of the 
coupling constant 2A /T C is close to the BCS value of 
3.53. 

We note that the obtained gap values A = 2.4(1) meV 
are in good agreement with the value of 2.5 meV de- 
rived by Kato et al. 37jj from scanning tunneling spec- 
troscopy of Fe1.05Seo.15Teo.85, by Homes et al. [38| from 
optical conductivity of FeSen4sTeo.55, and with 2.6 meV 
obtained by Biswas et al. [39( from /LtSR, and by Ben- 
dele et al. from magnetic penetration studies of 



FeSe .5Te .5. These studies and ARPES |4lJ indicate a 
multigap structure of FeSer-uTe^. A fit to the electronic 
specific heat using a two-gap model, however, did not 
give an essential improvement compared to the single- 
gap model and therefore we conclude that the specific 
heat is dominated by only one gap of 2.4 meV. 
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FIG. 11. (color online) Temperature dependences of the elec- 
tronic specific heat at different applied magnetic fields for 
oxygen-free sample F213 with bulk superconductivity. The 
arrow shows the direction of increasing field. 

Finally, in Fig. 11 we present the data for variation 
of the superconducting specific heat under applied mag- 
netic fields for one of the purest samples with high su- 
perconducting parameters. For better presentation the 
data are shown after subtracting the lattice contribution 
from the total specific heat. These dependences allow 
to get an estimate of the upper critical field H c2 from 
the temperature shift of the specific heat using the cri- 
terion of the minimum of the temperature derivative of 
C e in the transition region. The respective dependence 
H C 2 = f(H) obtained from these data is shown by open 
squares in Fig. 7. Interestingly, it exhibits a behavior 
resembling that of the H C 2 curve determined from the 
resistivity data for the field applied along the ab plane. 
The calculations by the WHH formula 2jJ gave a value 
of Hc2(0) ~ 1300 kOe which is by a factor of 2.5 larger 
than that obtained from the resistivity data. Similar re- 
sults were recently reported for Ba(K)Fe2As2 pnictides 



by Popovich et al. [42j], showing that H C 2{0) derived 



from the specific heat is by a factor of two higher than 
that calculated from the resistivity, which was attributed 
to flux-flow effect due to vortex motion. We also would 
like to mention that concerning H C 2 very similar data de- 
termined from specific heat and resistivity were recently 
reported by Serafin et al. (43[, which were interpreted 
as due to the presence of strong thermal fluctuations in 
FeSeo.5Teo.5- 
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TABLE III. Parameters determined from the heat capacity measurements 



Sample 


7r 

(mJ/mol K 2 ) 


a 

p 

(mJ/mol K 4 ) 


7n 

(mJ/mol K 2 ) 


A 

Ao 
(K) 


2A()/J C 


Br N5 


5.2 


0.75 


24 


26.6 


3.94 


F213 


0.82 


0.85 


25 


28.1 


3.86 


F216 step 1 


0.96 


0.94 


25 


25.9 


3.57 


F216 step 2 


19.3 


0.90 


23 







CONCLUDING REMARKS 

Our detailed studies of the structural, magnetic and 
thermodynamic properties of FeSeo.sTeo.5 reveal several 
important results: 

1. Preparation conditions have a substantial influence 
on the sample properties: The purity of the starting ma- 
terials and handling atmosphere are crucial for obtaining 
high-purity samples. Samples prepared from the puri- 
fied elements show superior properties compared to other 
cases. The oxygen-free samples prepared by fast cool- 
ing exhibit the lowest values of the susceptibility in the 
normal state, the highest transition temperature T° n of 
14.5 K, contain a volume fraction of the superconduct- 
ing phase up to 98%, and exhibit the most pronounced 
anomaly in the specific heat at T c . In the oxygen- free 
samples prepared by slow cooling the bulk superconduc- 
tivity is suppressed. They exhibit a paramagnetic tail 
in the susceptibility at low temperatures, a very small 
width of the hysteresis loop, and a strongly broadened 
resistive transition. 

2. The magnetic hysteresis measurements revealed 
high values of the critical current density j c of 8.6 x 
10 4 A/cm 2 for the purest samples which can be at- 
tributed to intrinsic inhomogeneity due to disorder at the 
anion sites. The oxidized samples show an increased j c 
up to 2.3 x 10 5 A/cm 2 due to additional pinning centers 
of Fe 3 4 . 

3. The upper critical field H c2 of ~ 500 kOe is esti- 
mated from the resistivity study in magnetic fields par- 
allel to the c-axis for both pure samples and samples 
containing oxide impurity. The anisotropy of the upper 
critical field -fH c2 = H^/H^ 2 reaches a value of about 6 
at T/T c = 0.996 and is the highest reported to date for 
these materials. 

4. The specific-heat measurements evidenced very low 
values of the residual Sommerfeld coefficient correspond- 
ing up to 96 % volume fraction of the superconducting 
phase which confirms the high quality of the oxygen-free 
samples. 

5. The temperature dependence of the electronic su- 
perconducting specific heat for samples with bulk super- 
conductivity can be reasonably described within a single- 
band BCS model with the temperature dependent gap 
A of value 27(1) K at T = K. The values of the cou- 



pling constant 2A /T C (see Tab. 3) are rather close to 
the BCS value of 3.53. 

6. The electronic specific heat of samples with sup- 
pressed bulk superconductivity shows a broad cusp-like 
anomaly which is ascribed to a splitting of the ground 
state of the Fe 2+ ions at the 2c sites. This contribution 
is fully suppressed in the ordered state in samples with 
bulk superconductivity. 

The authors thank Dana Vieweg and Nikola Pascher 
for experimental support. This research has been sup- 
ported by the DFG via SPP 1458 and Transregional Col- 
laborative Research Center TRR 80 (Augsburg - Mu- 
nich) . 



[1] Y. Kamihara, T. Watanabe, M. Hirano, H. Hosono, J. 

Am. Chem. Soc. 130, 3296 (2008). 
[2] M. Rotter, M. Tegel, D. Johrendt. Phys. Rev. Lett. 101, 

107006 (2008). 

[3] J.H. Tapp, Z. Tang, B. Lv, K. Sasmal, B. Lorenz, C.W. 
Chu, A.M. Guloy, Phys. Rev. B 78, 060505 (2008). 

[4] F.-C. Hsu, J.Y. Luo, K.W. Yeh, T.K. Chen, T.W. Huang, 
P.M. Wu, Y.-C. Lee, Y.L. Huang, Y.-Y. Chu, D.C. Yan, 
M.-K. Wu, Proc. Natl. Acad. Sci. U.S.A. 105, 14262 

(2008) . 

[5] J.W. Lynn, P. Dai, Physica C 469, 469 (2009). 

[6] S. Li, C. de la Cruz, Q. Huang, Y. Chen, J.W. Lynn, J. 

Hu, Y.-L. Huang, F.-C. Hsu, K.-W. Yeh, M.-K. Wu, and 

P. Dai, Phys. Rev. B 79, 054503 (2009). 
[7] S. Medvedev, T.M. McQueen, LA. Troyan, T. Palasyuk, 

M.I. Eremets, R.J. Cava, S. Naghavi, F. Casper, V. 

Ksenofontov, G. Wortmann, and C. Fclser, Nature Ma- 
terials 8, 630 (2009). 
[8] S. Margadonna, Y. Takabayashi, Y. Ohishi, Y. 

Mizuguchi, Y. Takano, T. Kagayama, T. Nakagawa, M. 

Takata, and K. Prassides, Phys. Rev. B 80, 064506 

(2009) . 

[9] T.M. McQueen, Q. Huang, V. Ksenofontov, C. Felser, Q. 
Xu, H. Zandbergen, Y.S. Hor, J. Alfred, A.J. Williams, 
D. Qu, J. Checkelsky, N.P. Ong, and R.J. Cava, Phys. 
Rev. B 79, 014522 (2009). 

[10] E. Pomjakushina, K. Conder, V. Pomjakushin, M. Ben- 
dele, and R. Khasanov, Phys. Rev. B 80, 024571 (2009). 

[11] M.H. Fang, H.M. Pham, B. Qian, T.J. Liu, E.K. Vehst- 
edt, Y. Liu, L. Spinu, and Z.Q. Mao, Phys. Rev. B 78, 
224503 (2008). 

[12] K.W. Yeh, T.-W. Huang, Y. Huang, T.-K. Chen, F.-C. 
Hsu, P.M. Wu, Y.-C. Lee, Y.-Y. Chu, C.-L. Chen, J.- 



11 



Y. Luo, D.-C. Yan, and M.-K. Wu, Europhys. Lett. 84, 
37002 (2008). 

[13] B.C. Sales, A.S. Sefat, M.A. McGuire, RY. Jin, D. Man- 
drus, and Y. Mozharivskyj, Phys. Rev. B 79, 094521 
(2009). 

[14] G.F.Chen, Z.G. Chen, J. Dong, W.Z. Hu, G. Li, X.D. 
Zhang, P. Zheng, J.L. Luo, and N.L. Wang, Phys. Rev. 
B 79, 140509 (R) (2009). 

[15] R. Khasanov, M. Bendele, A. Amato, P. Babkevich, A.T. 
Boothroyd, A. Cervelino, K. Conder, S.N. Gvasaliya, 
H. Keller, H.-H. Klaus, H. Luetkens, V. Pomjakushin, 
E. Pomjakushina, and B. Roessli, Phys. Rev. B 80, 
140511(R) (2009). 

[16] T. Taen, Y. Tsuchiya, Y. Nakajima, and T. Tamegai, 
Phys. Rev. B 80, 092502 (2009). 

[17] R. Hu, E.S. Bozin, J.B. Warren, and C. Petrovic, Phys. 
Rev. B 80, 214514 (2009). 

[18] J. Rodriguez-Carvajal, Physica B 192, 55 (1993). 

[19] E.J.W. Verwey, Nature (London) 144, 327 (1939). 

[20] The electron-spin resonance spectra of the samples han- 
dled in air show broad features characteristic of the Ver- 
wey transition observed in magnetite Fe304 which was 
investigated complimentarily. Below the Verwey transi- 
tion these features in the spectra fully disappear similar 
to the magnetite sample. The purest (oxygen-free) sam- 
ples do not show any noticeable features related to the 
Verwey transition. 

[21] T.J. Liu, X. Ke, B. Qian, J. Hu, D. Fobes, E.K. Vehstedt, 
H. Pham, J.H. Yang, M.H. Fang, L. Spinu, P. Schiffer, 
Y. Liu, and Z.Q. Mao, Phys. Rev. B 79, 174509 (2009). 

[22] B. Shen, P. Cheng, Z. Wang, L. Fang, C. Ren, L. Shan, 
and H.-H. Wen, Phys. Rev. B 81, 014503(2010). 

[23] CP. Bean, Phys. Rev. Lett. 8, 250 (1962). 

[24] CP. Bean, Rev. Mod. Phys. 36, 90 (1964). 

[25] R. Prozorov, N. Ni, M.A. Tanatar, V.G. Kogan, R.T. 
Gordon, C. Martin, E.C Blomberg, P. Prommapan, J.Q. 
Yan, S.L. Budko, and PC Canfield, Phys. Rev. B 78, 
224506 (2009). 

[26] J. Yang, M. Matsui, M. Kawa, H. Ohta, C. Michioka, C. 

Dong, H. Wang, H. Yuan, M. Fang, K. Yoshimura, arXiv: 

0911.4758 (unpublished). 
[27] N.R. Werthamer, E. Helfand, and PC Hohenberg, Phys. 

Rev. 147, 295 (1966). 
[28] B. Zeng, G. Mu, H.Q. Luo, T. Xiang, H. Yang, L. 

Shan, C. Ren, I. I. Mazin, P.C Dai & H.-H. Wen, 

larXiv:1004.2236l (unpublished). 
[29] V. Ksenofontov, G. Wortmann, A.I. Chumakov, T. Gasi, 

S. Medvedev, T.M. McQueen, R.J. Cava, and C. Felser, 

Phys. Rev. B 81, 184510 (2010). 
[30] D. Phelan, J.N. Millican, E.L. Thomas, J.B. Leao, Y. 

Qiu, R. Paul, Phys. Rev. B 79, 014519 (2009). 



[31] In several earlier papers on FeSe(Te) a fit by an odd- 
power polynomial to the experimental data taken just 
above T c was used to separate the electronic and lattice 
specific heat. Utilizing a similar fitting procedure for the 
temperature range 15 - 21 K we obtained for the nor- 
mal electronic coefficient value 7„ = 90 mj/mol K 2 and 
for the prefactor /3 = 0.3 mj/mol K 4 corresponding to 
a Debye temperature Qd ~ 235 K. However, the sim- 
ple Debye approximation is known to work well only for 
temperatures below Od/50 = 4.7 K [32|, which is much 
lower than the temperature range of fitting. Therefore 
these values of 7„ and 9d, are strongly overestimated. 
The failure of the harmonic-lattice approximation above 
T c was demonstrated recently also for the related super- 
conducting BaFe(Co)2As2 pnictides [33T] . 

[32] C. Kittel, Introduction to Solid State Physics, John Wi- 
ley and Sons, Inc. New York, London, Sydney, Toronto, 
1971. 

[33] F. Hardy, T. Wolf, R.A. Fisher, R. Eder, P. Schweiss, 

P. Adehnann, H.v. Loehneysen, and C. Meingast, Phys. 

Rev. B 81, 060501(R) (2010). 
[34] F. Bouquet, Y. Wang, R.A. Fisher, D.G. Hinks, J.D. 

Jorgensen, A. Junod, and N.E. Phillips, Europhys. Lett. 

56, 856 (2001). 

[35] H. Padamsee, J.E. Neighbor, and C.A. Shiffman, J. Low 

Temp. Phys. 12, 387 (1973). 
[36] Ch. Kant, J. Deisenhofer, A. Giinther, F. Schrettle, A. 

Loidl, M. Rotter, and D. Johrendt, Phys. Rev. B 81, 

014529 (2010). 

[37] T. Kato, Y. Mizuguchi, H. Nakamura, T. Machida, H. 
Sakata, and Y. Takano, Phys. Rev. B 80, 180507 (R) 

(2009) . 

[38] C.C. Homes, A. Akrap, J.S. Wen, Z.J. Xu, Z.W. Lin, Q. 

Li, and CD. Gu, Phys. Rev. B 81, 180508 (2010). 
[39] PK. Biswas, G. Balakrishnan, D.M. Paul, C.V. Tomy, 

M.R. Lees, and A.D. Hillier, Phys. Rev. B 81, 092510 

(2010) . 

[40] M. Bendele, S. Weyeneth, R. Puzniak, A. Maisuradze, E. 

Pomjakushina, K. Conder, V. Pomjakushin, H. Luetkens, 

S. Katrych, A. Wisniewski, R. Khasanov, and H. Keller, 

Phys. Rev. B 81, 224520 (2010). 
[41] F. Chen, B. Zhou, Y. Zhang, J. Wei, H.-W. Ou, J.-F. 

Zhao, C. He, Q.-Q. Ge, M. Arita, K. Shimada, H. Na- 

matame, M. Taniguchi, Z.-Y. Lu, J. Hu, X.-Y. Cui, D.L. 

Feng, Phys. Rev. B 81, 014526 (2010). 
[42] P. Popovich, A.V. Boris, O.V. Dolgov, A. A. Golubov, 

D.L. Sun, C.T. Lin, R.K. Kremer, and B. Keimer, Phys. 

Rev. Lett. 105, 027003 (2010). 
[43] A. Serafin, A. I. Coldea, AY. Ganin, M.J. Rosseinsky, 

K. Prassides, D. Vignolles, and A. Carrington, arXiv: 

1007.3914 (unpublished). 



